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ABSTRACT 

We present a study of the Angular Resolution of the AGILE gamma-ray imaging 
deteetor (GRID) that is operational in spaee sinee April 2007. The AGILE instrument 
is made of an array of 12 planes eaeh equipped with a Tungsten eonverter and Silieon 
mierostrip deteetors and is sensitive in the energy range 50 MeV - 10 GeV. Among 
the spaee instruments devoted to gamma-ray astrophysies, AGILE uniquely exploits 
an analog readout system with dedieated eleetronies eoupled with Silieon deteetors. 

We show the results of Monte Carlo simulations earried out to reproduee the gamma- 
ray deteetion by the GRID, and we eompare them to in-flight data. We use the Crab 
(pulsar -I- Nebula) system for diseussion of real data performanee, sinee its E~^ energy 
speetrum is representative of the majority of gamma-ray sourees. For Crab-like spee- 
trum sourees, the GRID angular resolution (FWHM of ~ 4° at 100 MeV; ~ 0.8° at 1 
GeV; ~ 0.9° integrating the full energy band from 100 MeV to tens of GeV) is stable 
aeross a large field of view, being eharaeterized by a flat response up to 30° off-axis. 

A eomparison of the angular resolution obtained by the two operational gamma-ray 
instruments, AG/LE-GR/D and Eermi-LAT, is interesting in view of future gamma-ray 
missions, that are eurrently under study. The two instruments exploit different deteetor 
eonfigurations affeeting the angular resolution: the former being optimized in the read¬ 
out and traek reeonstruetion espeeially in the low-energy band, the latter in terms of 
eonverter thiekness and power eonsumption. We show that, despite these differenees, 
the angular resolution of both instruments is very similar between 100 MeV and a few 
GeV. 

Subject headings: gamma rays: instruments, gamma rays: sourees 


1. Introduction 


Gamma-ray astrophysies in spaee enormously advaneed sinee the first deteetion of photons 
above 100 MeV (iKraushaar et al. 19721) . The progression of spaee missions and instruments over 
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the last decades (OSO-3, SAS-2, COS-B, EGRET on board of the Compton Gamma-Ray Observa¬ 
tory) led to improvements of the overall detector performance, in terms of both angular resolution 
and sensitivity. In addition to this, the progressively wider Eield-of-View (EoV) allowed a contin¬ 
uous monitoring of the variable gamma-ray sk y. Eollowing the early cosmic gamma-ray detection 
by the OSO-3 satellite (IKraushaar et al. 19721) . the fir st gamma-ray teles cope SAS-2 launched in 
1972 reached an angular resolution of a few degrees (lEichtel et al. 1975h . This mission was fol- 
lowed by the European mission COS-B, launched in 1975 August 8 and operational for 7 years 
dMayer-Hasselwander et al. 19791;ISwanenburg et al. 19811) . The Compton Gamma-Ray Observa¬ 
tory (CGRO), active between 1991-2000, provided the first complete investigation of the gamma- 
ray sky. In particular, CGRO hosted the Energet ic Gamma-Ray Experiment Telescope (EGRET), 


operating in the energy range 30 MeV - 30 GeV (lEichtel & Trombka 1997l:lThompson et al. 1993 


Thompson et al. 1998h . Fre-AGILEfPre-Fermi space instruments were mostly based on spark cham¬ 
ber technology. 

Eurther improvements to the overall performance of gamma-ray detectors in space became 
possible with the advent of solid-state Silicon detector technology. The scientific objectives of the 
new generation of instruments required the following enhancements: (1) improving the gamma-ray 
angular resolution near 100 MeV by at least a factor of 2-3 compared to EGRET; (2) obtaining the 
largest possible EoV at 100 MeV reaching 2.5-3 sr ; (3) increasing flux sensitivity near 100 MeV. 
The current generation of gamma-ray space instruments, AG/LE-GR/Z) (Gamma-Ray Imaging De¬ 
tector) and Fermi-LAT (Earge Area Telescope), launched in April 2007 and June 2008 respectively, 
were designed to achieve these objectives. Both instruments make use of tungsten-silicon detec¬ 
tors for the conversion and detection of gamma-ray photons , having a common o verlapping energy 
range in the ~ 50 MeV - 10 GeV band ( Tavani et al. 2009 : Atwood et al. 2009 ). The gamma-ray 
detector is structured to form a “Silicon Tracker” made of several trays, each containing a tung¬ 
sten layer (used as converter of the incident y-ray photon into an e'^/e“ pair) and two silicon strip 
layers (used to track the the e'^/e“ pair path across the instrument, through their (x,y) projections 
along each tray). The detected tracks are identified and fitted by a Kalman filter that allows for 
the reconstruction of the original direction and energy of the incident photon. The main source 
of contamination in the detection of y-ray photons is then by charged particle tracks, which can 
be confused with pairs. Eor this reason, background rejection filters are used both on-board 
and on-ground, to process the data and obtain a final discrimination and classification of the events. 


An important parameter in the assessment of the overall instrument scientific performance 
is the resulting angular resolution, i.e. the minimum distance at which two close sources are 
distinguished as separated. The Point Spread Eunction (PSE) describes the response of an imaging 
instrument to a point source and the angular resolution is usually described by either the Eull Width 
at Half Maximum (EWHM) of the PSE radial profile or the 68% containment radius (CR68%) of 
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the PSF. Key parameters of the instrument eonfiguration affeeting the angular resolution of y-ray 
solid-state silieon deteetors are: 


• the size of the silieon strips (the piteh) and the silieon deteetor readout system (e.g., either 
digital or analog) resulting in different effeetive spatial resolution of the partiele trajeetory 
due to eharge coupling between adjacent silicon strips; 

• the distance between consecutive silicon planes that, combined with the spatial resolution, 
defines the limiting angular resolution; 

• the thickness of tungsten layers promoting gamma-ray photon conversion, and at the same 
time inducing multiple scattering that leads to a degradation of the charged particle tracks; 

• the reconstruction and event classification algorithms. 


Table A1 in the Appendix summarizes the main characteristics of the AGILE and Fermi 
gamma-ray detectors, whose configurations are different in several ways: the former being op¬ 
timized in the readout and track reconstruction especially in the low energy band, the latter in 
terms of converter thickness, geometrical area, and power consumption. An important difference 
between the two instruments, that can be crucial for the scientific performance of the instrument, 
is the readout system (analog for iht AGILE-GRID and digital for Fermi-LAT). 


In the following sections, we focus on the characterization of the AGILE-GRID angular reso¬ 
lution both from simulations and from in-flight data. In section [2] we define the parameters that are 
used across the paper to describe the instrument angular resolution; section [3] describes the Monte 
Carlo (MC) simulations setup and data processing pipeline; section |4] shows the results from the 
MC simulations; section [5] compares the results of the simulation with the in-flight data angular 
resolution. This paper complements the work by dChen et al. 2013h : here we focus on simulations 
characterizing the overall performance of the GRID instrument and in a special study of Crab-like 
sources, reproducing the behaviour of the majority of detected cosmic gamma-ray sources; we also 
present a direct comparison of AGILE and Fermi in-flight data for the CrabQ, and discuss the results 
concerning the angular resolution in terms of the different instrument configurations (section [^, 
which can be crucial in the study of future missions. 


'The case of the Vela pulsar has been addressed by dChen et al. 2013h and (lAckermann et al. 2013h . 
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2. The Angular Resolution 

For an imaging telescope, the response in terms of reconstruncted positions of a set of photons 
from a point-like source in the sky defines its Point Spread Function (PSF). Since this function can 
be considered azymuthally symmetric in y-ray telescopes for incidence angle within 30° as used in 
this paper, the effect of the dispersion can be described as a function of one parameter, the angular 
distance a between reconstructed and nominal direction. Defining the PSF radial profile, PCa), 
as the probability distribution per steradian of measuring an incoming photon at a given angular 
distance a from its true direction, we have therefore: 


PSF(a)da=27Tsin(a)P(a)da. 

In the following, we adopt two parameters w idely used in the litera ture to assess the angular 
resolution in terms of the “width” of the PSF (e.g. lThompson et al. 1993 ): 

• the Full Width at Half Maximum (FWHM) of the probability distribution per steradian, P(a) 

• the 68% Containment Radius (CR68%) of the probability distribution, PSF(Qr). 

The CR68% is strongly related to the source image compactness, taking into account the whole 
contribution of source profile and being more affected by possible extended tails; while the FWHM 
is mainly determined by the central core emission of the source. 


3. Setup and data processing of the AGILE simulated data 


In order to fully understand how the instrument configuration and the data analysis pipeline 
affect the final angular resolution of a gamma-ray telescope, we carried out dedicated MC sim¬ 
ulations of the AGILE-GRID analyzing simulated data with the same pipeline used to process 
in-flight data. The simulations were carried out using an available s imulation tool (“GAMS”, 
GEANT AGILE MC Simulator; Cocco et al. 2002 : Congo et al. 2002 ). implemented during the 
development phase of the AGILE mission. The tool makes use of the GEANT3 environmenu and 
takes into account all the main components of the AGILE instrument configuration: the space¬ 
craft (bus) MITA and the AGILE payload, consisting of the Csl Mini-Calorimeter, the silicon- 
tungsten Tracker, the Anti-Coincidence system, the X-ray detector (Super-AGILE), the thermal 


^ http://wwwasd.web.cern.ch/wwwasd/geant/index.html 
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shield, the mechanical structure and the lateral electronics boards. The tool allows to simulate a 
parallel front of given direction for both charged particles and/or photons and its interactions across 
the instruments. The spectral energy distribution of the front can be monochromatic or with any 
given law. The simulated tracks are t hen processed by the DHSIM, the Data Handling SIMulator 
(lArgan et al. 20041: lArgan et al. 20081) . which implements the on-board algorithms for a first track 
reconstruction, event classification and background rejection (IGiuliani et al. 20061) . The DHSIM 
output is a data file in the same format as the AGILE in-flight telemetry and can be processed using 
the same pipeline as for real data analysis. Further background rejection techniques are applied to 
the data on ground, producing the final classification of the events. The currently used filter for 
the scientific analysis of the AGILE-GRID data (“FM3.119”) was used in this paper also for the 
processing of the simulated data. 

Each simulation contains 2x10^ photons that cross the tracker and the data are analyzed with the 
same pipeline as for real data. Typical efficiency after background rejection within 30° off-axis 
angle is ~ 30%, leaving therefore a set of ~ 5 x 10^ events classified as photons and used for our 
calculations in the paper. 


4. The Angular Resolution for simulated AGILE-GRID data 

With the aim of investigating the dependence of the telescope angular resolution upon the 
energy of incoming photons and their incident direction in the FoV, we carried out simulations of 
parallel fronts of photons by varying those parameters. 

We describe the incident direction as the composition of a zenith angle 6 (i.e. the angular distance 
of the incoming photon direction from the vertical axis of the Tracker, also named off-axis angle) 
and an azimuthal angle (p (be., the angular distance between a given axis in the instrument plane 
and the projection of the incoming direction in the plane). Since the detector response is mostly 
azimuthally symmetric, we expect the main dependence to be upon 6. Simulations were carried 
out for the following values of these angles: 6 = (1°, 30°, 50°)§; (/> = (0°, 45°). Two different sets of 
simulations were performed: 

• parallel fronts of monochromatic photons of energies 50, 100, 200, 400, 1000, 5000 MeV; 

• parallel fronts of photons with a Crab-like photon spectrum of the type \ in the energy 
range 30 MeV - 50 GeV. 


^For the on-axis case we use a value of 6 = 1°, since the case of 0 = 0° induces a singularity in the software that 
generates the parallel front of photons. For our purposes 0=1° well resembles the on-axis case. 
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The simulated data are analysed with the same pipeline as for the real data. In partieular data 
were analysed using the most reeent pipeline for the data proeessing (“BUILD21”) and instrument 
response funetions (“10023”). 
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■■ 

140 - 300 

1.1 

1.2 
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0.1 
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0.2 
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Table 1: Results from the simulations of monochromatic photon parallel beams. Half Width at Half Maxi¬ 
mum of the PSF radial profile in degrees. Events for the different energy channels are selected on the basis 
of the reconstructed energy, as for real data. 
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Table 2: Results from the simulations of monochromatic photon parallel beams. 68% containment radius 
of the PSF in degrees. Events for the different energy channels are selected on the basis of the reconstructed 
energy, as for real data. 
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4.1. Monochromatic photons 

Monochromatic simulated data can be used for an ideal characterization of the instrument 
response and were simulated for different energies and different incident directions as previously 
described. 

According to the definitions given in section [21 Table [H and [2] show the values obtained for 
the HWH1V0 (calculated from the PSF radial profile) and for the CR68% (calculated from the PSF) 
for different energies and different off-axis angles. We use both parameters for symplifing the 
comparison with previous papers, related both to gamma rays or to other wavelengths. In the rest 
of the paper we prefer adopting the FWHM to describe the angular resolution of the instruments 
in order to be more compliant with the definition in the multi-frequency domain. 

Note that for this analysis we do not fit the PSF with any specific function, we just obtain the 
above parameters from the raw distribution of the reconstructed directions (see e.g. Fig. [fright 
panel). The source radial profile is given by the average source counts within circular crowns of in¬ 
creasing radii, centered at the source centroid. The instrumental background counts are subtracted 
by evaluating their contribution in a circular crown at distance large enough to avoid contamina¬ 
tion by the source tails (see e.g. Fig. [S]). In Tables [Hand [2] we also report the intrinsic error in the 
estimate, given by the bin size of the profile, which has a minimum value of 0.1°. This value is 
the minimum bin size allowed for the generation of typical count maps in the AGILE real data for 
high significance sources. At the same time this value allowed us to have a statistical significance 
well above 5-cr per bin in the radial profiles used to calculate the FWHM with the adopted number 
of simulated input events. 

Table[I]shows that iht AGILE angular resolution, as inferred from the FWHM, is stable across 
the instrument field of view up to 30° of off-axis angle. Table [2| shows that, on the contrary, 
the CR68% increases slightly with off-axis angle below 1000 MeV, since this parameter is much 
more affected by tails in the radial profile than the FWHM. This is expected, since at 30° off-axis 
angle tracks pass through much more material compared to the on-axis case: this causes a higher 
dispersion in the reconstructed direction due to the effect of multiple scattering, more pronounced 
below 1000 MeV (see Sec. [6|). 


"^Here we show the HWHM in place of the FWHM used in the figures in order to help for a direct comparison with 
the CR68% and in the identification of the width of the radial profile (see e.g. Fig[2l right panel). 
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Fig. 1.— Left panel: AGILE-GRID 68% containment radius versus photon energy for simulated 
monoehromatie photons of different ineident angles. Right panel: angular resolution (FWHM) 
versus photon energy, shown for eomparison to the CR68. 


4.2. Simulations of Crab-like sources 

Simulations of Crab-like speetrum sourees were earned out as previously deseribed and, as 
an example. Fig. |2] shows a eount map obtained from the MC simulated data for the 100 - 400 
MeV energy band at the (30°, 0°) ineident direetion (left panel) and the average radial profile of 
the souree (right panel). The map was obtained by the simulated event list using the same pipeline 
as that for produeing real data maps. 

Table [3] shows the results of the whole analysis for the Crab simulations at varying photon energy 
and ineident direetion. The instrument response is stable within 30° aeross the FoV within the 
errors and the overall response is dominated by the zenith angle (6), sinee the dependency upon 
azymuthal angle is minimal, if any. Fig. [2] shows the angular resolution as a funetion of photon 
energy for on-axis ineident direetions. 

We also report in Table [3] the HWHM for the typieal broad bands used in the AGILE-GRID data 
analysis, i.e. the 100 - 400 MeV, the 400 - 1000 MeV ehannels and the 100 MeV - 50 GeV full 
band, whieh gives on average an angular resolution of 1.8° + 0.2°, 0.8° + 0.2° and 0.9° + 0.2° for 
off-axis angles within 30° respeetively. These values will be compared to the ones obtained with 
the in-flight data analysis. 

It is interesting to compare the results of the Crab simulations to the monoehromatie photon 
beams (see Tabled and [U): the values of the angular resolution per energy ehannels agree within 
the errors, showing that the eontribution of possible energy ehannel eross talk is negligible. 
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Fig. 2.— Left panel. Count map of the simulated Crab system in the 100 - 400 MeV energy band. 
Pixel size is 0.1°. Right panel Average eount radial profile. The eounts are normalized to the 
maximum value of the profile in order to reaeh 1 at 0°. The error bars are the poissonian values. 
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Table 3: HWHM of AGILE-GRID PSF radial profile for different energy ehannels and off-axis 
angles for the Crab simulations. 


5. In-flight data analysis 

The whole simulation setup and the above results for the angular resolution ean be validated by 
a eomparison with the values obtained by in-flight data. The AGILE satellite has been operational 
sinee July 2007 and the angular resolution of the GRID has been studied both on-ground (Cattaneo 
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Fig. 3.— AGILE-GRID angular resolution versus photon energy for the simulation of an on-axis 
Crab-like souree. 

et ah, 2012; Cattaneo & Rappoldi, 2013) and in-flight (Chen at el., 2013). Here however we foeus 
in partieular on the case of the Crab source (pulsar -l- Nebula), whose spectral energy distribution 
in the gamma-ray energy band is dominated by the pulsar and is described by a power law with a 
spectral index of ~2.1. As previously mentioned, this spectral energy distribution is representative 
of the majority of the gamma-ray sources and we therefore think that it is an ideal test case to assess 
the instrument response, also in view of future gamma-ray missions under study (see next section). 
We selected AG/LE data from a long-term pointing of the region (11 days, 2007-09-23X12:01:05 - 
2007-10-04X12:01:05), during which the source was located at an off-axis angle < 30°. Due to the 
resulting overall exposure of the region, we focus hereafter on two wide energy bands in order to 
obtain a good photon statistics for a robust assessment of the angular resolution: the 100 - 400 MeV 
and 400 MeV - 1 GeV energy bands. Fig. [5] and 6 show the Crab intensity map for the two energy 
bands respectively and the average count radial profile for the source. Xhe FWHM associated to 
these profiles turns out to be 2.5° + 0.5° in the 100 - 400 MeV energy range, and 1.2° + 0.5° for the 
400 MeV - 1 GeV band. Xhe good agreement with the results from the simulations (see XableO 
validate the whole setup and data analysis carried out for the simulations. 
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Fig. 4.— Left panel: AGILE-GRID count map of the Crab system (pulsar + Nebula) for the period 
2007-09-23X12:01:05 - in 2007-10-04X12:01:05 in the 100 - 400 MeV energy band; the white 
eirele shows the maximum aperture (6°) used to produee the radial profile of the souree shown 
in the right panel; green eireles show the erown region used to estimate the baekground emission 
(7.5° - 8.5° radius). Right Panel: Average eount radial profile: blue data points are the in-flight 
data, while the gray shaded region shows the profile for the simulated erab system with errors. Xhe 
eounts are normalized to the maximum value of the profile in order to reaeh 1 at 0°. Xhe souree in 
the up left eorner is Geminga. 


6. AGILE and Fermi: comparison of in-flight data for the Crah system 


As already mentioned in seetionfH although implementing the same silieon-tungsten deteetor 
eoneept, the instrument eonfigurations for the AGILE-GRID and Eermi-LAT deteetors are quite 
different in terms of eonverter thiekness, distanee between trays and spaeing between adjaeent 
silieon strips (see Xable Al). In the ease of the AGILE-GRID, the Silieon Xraeker is eomposed 
of 10 tungsten eonverter planes of homogeneous thiekness (0.07 radiation length eaeh) plus 2 
additional planes without eonverter (for a total of 12 planes); the overall GRID radiation length 
is ~ IXo dBarbiellini etal. 2002 : Prest et al. 2003 ). Xhe Eermi-LAT has 12 tungsten planes of 
~ 0.03 radiation length thiekness (“LAX-front”), and an additional s et of 4 tungsten plan es of 0.18 
radiation length (“LAX-baek”) plus two planes without eonverter (lAtwood et al. 20091) . In both 
instruments, eaeh tungsten plane is interleaved with 2 layers of silieon strip deteetors whieh are 
sensitive to eharged partieles and are used for the x,y positioning of the e'^/e“ pair traek resulting 
from the pair eonversion. Xhe ratio of strip piteh to vertieal spaeing between traeker planes is 0.007 
for both Eermi-LAT and AGILE-GRID. 
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Fig. 5.— Left panel: AGILE-GRID count map of the Crab system (pulsar + Nebula) in the 400 - 
1000 MeV energy band. Right Panel: Average eount radial profile: blue data points are in-flight 
data, while the gray shaded region shows the profile of the simulated erab system with errors. The 
eounts are normalized to the maximum value of the profile in order to reach 1 at 0°. 


An important difference between the two instruments, is the Silieon traeker readout system, 
analog for AGILE-GRID and digital for Eermi-LAT. Barbiellini et al. (2002) and Prest et al. (2003) 
deseribe extensively the GRID dedieated front-end eleetronies, eharaeterizing the effective spatial 
resolution for partiele ineidenee between 0° and 30° and eomparing the performanees of digital to 
analog readout. Although the GRID analog readout is structured to read only odd-numbered strips 
with no signal piek up at even-numbered strips (“floating strip readout”), the eapaeitive eoupling 
between adjaeent strips allows to obtain a complete sampling of the partiele hit positions, discrim¬ 
inating between hits involving direetly read and not-read mierostrips. Typically, each particle hit is 
charaeterized by a signal spread out over several adjacent readout-strips (a “cluster”). GRID events 
produce energy deposition histograms per eluster whieh are very well eharacterized and typieally 
result in 2-3 triggered strips depending on off-axis angle (Barbiellini et al. 2002; Fedel et al. 2000). 
Note that also for on-axis ineidenee a eluster is eomposed of a minimum of 2 strips. The effeetive 
spatial resolution obtained by the GRID for partiele ineidenee between 0° and 30° is 6, ~ 40 urn , 
substantially less than the Silieon mierostrip size of 121 pm (see Fig. 15 in [Barbiellini et al. 2002h . 
The advantage of the an alog vs. digital readout for the resulting spatial resolution of silicon strip 
detectors is discussed in (IBarhiefiini et al. 2002h for the AGILE traeker eonfiguration, showing an 
improvement of a faetor of 2 (i.e. 6s ~ 40 pm vs. 6s ~ 80 pm). 

The Eermi-LAT system is based on a digital readout. The EAT Traeker is non-homogeneous 
and is eharaeterized by two different values for the angular resolution, the one for the “EAT-front” 
and the of the “EAT-back” (lAekermann et al. 2013h . In the following, we use only EAT-front data 
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Fig. 6.— Comparison of the expected limiting values for the angular res¬ 
olution (dot-dashed cyan line) to the measured ones (red points). Left 

panel: AGILE-GRID simulated Crab data; right panel: Fermi-LAT PSF68 

(http://www. slac. Stanford. edu/exp/glast/groups/canda/lat_Perforinance .htm). 
Top panels. The overall limiting angular resolution is plotted in cyan and is obtained by the 
combination of: 1. the ratio of the spatial resolution to the distance of consecutive trays (dashed 
blue line); 2. the multiple scattering limit (red dot-dashed line). Bottom panels. Percentage 
Deviation of the measured angular resolution from the overall expected limiting value. 

for comparison with the AGILE-GRID, proving the best angular resolution for Fermi. Note also 
that Fermi usually operates in scanning mode, with sources observed most of the time at large off- 
axis angles, with a consequent degradation of the angular resolution for a given source, compared 
to the on-axis performance. 

In order to compare the angular resolution of the two instruments at similar conditions, we 
selected Fermi data of the Crab system in pointing mode. We analyzed therefore a comparable set 
of data, with an exposure ensuring a similar photon statistics and with both gamma-ray detectors in 
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pointing mode. The AG/LE data in pointing modqj used for this comparison are the ones described 
in section [51 The Fermi data consist of 4 days of integration in the period 2012-07-04T23:24:44 
- 2012-07-08T10:44:43 during which the satellite was stably pointing at the Crab following a 
gamma-ray flaring episode, with the source located at 10° of off-axis angle in the FoV. We used 
the best available quality cut for the Fermi data, (the“ULTRA-CLEAN” event class, LAT-front 
photons, P7REP) in order to optimize the angular resolution of the data. We focus on the analysis of 
two representative energy bands for AGILE in order to have a robust statistics for both instruments: 
100 - 400 MeV, and 400 - 1000 MeV. 

Pig. |7] shows the intensity maps for the two gamma-ray telescopes in the two energy bands 
and Pig. [8] shows the corresponding radial profiles. The resulting angular resolution deduced from 
the in-flight data are the same, even though the AGILE data show a more pronounced non-gaussian 
tail, due to the converter thickness: the FWHM are 2.5° + 0.5° and 1.2° + 0.5° in the energy range 
100 - 400 MeV and 400 MeV - 1 GeV respectively. 


As already mentioned, several factors determine the instrument angular resolution as a func¬ 
tion of gamma-ray photon energy, including: (1) multiple scattering, (2) the effective spatial reso¬ 
lution, (3) the photon energy reconstruction, (4) specific particle track reconstruction algorithms, 
(5) quality cuts. Pig. [6] shows the relative contibution of these factors as a function of energy 
for both trackers, AGILE and Fermi. Regarding the multiple scattering, we assume an effective 
radiation length per tray of 0.085 and 0.04 for AGILE-GRID and Fermi-LAT respectively, taking 
into account the contributions of the Tungsten converter and the supporting materia]^. Top panels 
show that the multiple scattering effect dominates the overall angular resolution up 700 MeV and 
350 MeV respectively for the two instruments. Bottom panels show the percentage deviation of 
the measured angular resolution to the overall expected limiting value. The amplitude of the devia¬ 
tion quantifies the reconstruction accuracy, that is optimized in different energy ranges for the two 
instruments: the accuracy is below 40% up to 400 MeV in the case of AGILE and above 2 GeV 
in the case of Fermi. Besides the different tracker configuration, crucial features distinguishing 
the two instruments and affecting the final reconstruction accuracy are the effective spatial resolu¬ 
tion and the quality of the reconstruction algorithmjzi Therefore w e believe that in view of future 
gamma-ray missions (e.g. lMorselli et al. 2013 : Galper et al. 2013 ) the optimization of the angular 


^An analysis of AGILE-GRID data in “spinning mode” will be discussed in a forthcoming paper by Lucarelli et 
al., 2015. 

®In the case of AGILE the effective radiation length due to the supp orting material for e ach tray is ~ 0.015 V 
from laboratory measurements; since Fermi has a similar tray structure (lAtwood et al. 2009h we assume a similar 
contribution of ~ 0.01 V- 


^Note that ongoing developments for the track reconstruction algorithms and events classification for both AGILE 
(updated FM hlter) and Fermi-LAT (PASS8) may lead to further improvements of the current angular resolution. 
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resolution should not only rely on a redueed eonverter thiekness, but also on the implementation 
of an analog readout system if the power budget allows it, together with an optimization of the 
reeonstruction algorithms. 



Fig. 7.— AGILE-GRID (left panels) and Fermi-LAT (right panels) in-flight count maps for the 
Crab region in the 100 - 400 MeV and 400 - 1000 MeV energy bands respectively. Pixel size is 
0.25°. Note that the color intensity scale is set so that each map has the maximum value of the 
scale corresponding to the Crab peak counts. 


7. Conclusions 

The current generation of gamma-ray space instruments is based on Silicon detector tech¬ 
nology and associated electronics. Both AGILE-GRID and Fermi-LAT show a quite stable per¬ 
formance in orbit which is the basis for prolonged operations {AGILE is at its 9th year of life in 












- 16- 


100-400 MeV 400-1000 MeV 




Radial Aperture (deg) Radial Aperture (deg) 

Fig. 8.— Average eount radial profiles for cireular aperture of inereasing radii at steps of 0.25° of 
the Crab (pulsar + Nebula); in-flight data for the AGILE-GRID (blue data points) and Fermi-LAT 
(red data points). Left panel: 100 - 400 MeV energy range; right panel: 400 - 1000 MeV energy 
range. 


spaee, Fermi at its 8th year). Compared to the previous generation, the instrument performanee 
improvements both in terms of sensitivity and angular resolution is well established. 

In this paper, we summarized the main results concerning the GRID angular resolution, a 
crucial feature of the AGILE instrument, at the basis of the scientific performance of a gamma-ray 
detector together with its FoV and background rejection capabilities at energies below 400 MeV, 
allowing for the best exploitation of the instrument configuration. We showed that the AGILE- 
GRID angular resolution is optimized given the overall characteristics of the detector and allows 
for state-of-the-art mapping of Galactic and extragalactic regions/sources: the FWHM for off-axis 
angles in the range 0° - 30° is ~ 4° at 100 MeV, and ~ 0.8° at 1 GeV. The angular resolution is 
quite uniform in the FoV up to 30° off-axis. The FWHM obtained from in-flight data in pointing 
mode is 2.5° in the range 100 - 400 MeV, and ~ 1.2° in the 400 MeV - 1 GeV. 


Although the AGILE-GRID multiple scattering due to the heavy converter is relatively high, 
we proved that a crucial role is played by the optimization of the readout system (analog) of the 
Silicon Tracker and of the particle track reconstruction algorithm^. 


is interesting to note that the AGILE Tracker configuration is quite similar to the basic element of the gamma-ray 
instrument currently under study for the GAMMA-400 mission (Galper etal. 2013). The analog readout of a Silicon 
Tracker with AGILE-like characteristics is required to optimize the angular resolution with a thick converter (that in 
the case of GAMMA-400 is currently designed to be ~ 0.08 V per plane). 
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The GRID angular resolution as a function of gamma-ray energies is shown in Fig. [T] as 
resulting from simulations. These values are in good agreement with the ones deduced from in¬ 
flight data as demonstrated in this paper for Crab-like sources. 

Furthermore, by a direct comparison of in-flight data of the Crab system, we find that de¬ 
spite the differences in structure, geometry, readout system and overall size, the AGILE-GRID and 
Fermi-LAT front show similar angular resolutions at energies between 100 MeV and 1 GeV, due to 
different optimizations of the readout system and reconstruction algorithms. 


Acknowledgements 

We ackowledge several discussions with our colleagues of the Fermi Team. The AGILE mis¬ 
sion is funded by the Italian Space Institute (ASI) with scientific and programmatic participation 
by the Italian Istitute of Astrophysics (INAF) and the Italian Insitute of Nuclear Physics (INFN). 
Our research is partially supported by the ASI grants 1/042/10/0,1/028/12/0 and 1/028/12/02. We 
would like to thank the referee for the careful review and for providing valuable comments that 
helped in improving the contents of this paper. 


REFERENCES 

Ackermann et ah, 2013, ApJ, 765, 19 

Argan, A. et ah, 2004, Proceedings lEEE-NSS, vol. 1,371. 

Argan A., Tavani M., Trois A., et ah, 2008, NSS. IEEE (IEEE), 774 
Atwood W.B. and the Eermi EAT collaboration, 2009, Apl, 697, 1071 
Barbiellini G., Eedel G., Eiello E. et al., 2002, NIM A, 490, 146 

Cattaneo P.W., Argan A., Boffelli A., Buonomo B., Chen A.W., D’Ammando F, Foggetta E., 
Eroysland T. et al., 2012, NIM A, 674, 55 

Cattaneo P.W., Rappoldi P, 2013, NIM B, 239, 241 

Chen A.W. et al., 2013, A&A, 558, 37 

Cocco V., Eongo E, Tavani M., 2002, NIM A, 486,623 

Eedel G. et al., 2000, SPIE, Vol. 4140 


18- 


Feroci M. et al, 2007, NIM A, 581, 728. 

Fichtel C.E., Hartman, R.C., KnifFen, D.A., Thompson, D.J., Bignami, G.F., Ogelman, H.B., Ozel, 
M.E., and Turner, T., 1975, ApJ, 198, 163. 

Eichtel C.E. & Trombka J.I., 1997, Gamma-Ray Astrophysics, NASA Reference Publication n. 
1386, September 1997. 

Galper A.M. et al, 2013, Proceedings of the International Cosmic-Ray Conference 2013, Brazil, 
Rio de Janeiro, http://arxiv.org/ftp/arxiv/papers/1306/1306.6175.pdf 

Giuliani A., 2003, The Kalman Algorithm for the AGILE Mission, Ph.D. Thesis, University of 
Milan 

Giuliani A., Cocco V., Mereghetti S., Pittori C., Tavani M., 2006, NIM A, 568, 692 

Hartman R.C. et al., 1999, ApJS, 123, 279 

Kraushaar W.E. et al., 1972, ApJ, 177, 341 

Eabanti C. et al., 2006, proc SPIE, 6266, 62663 

Eongo E, Cocco V., Tavani M., 2002, NIM A, 486, 610-622 

Eucarelli E. et al., 2015, in preparation 

Mayer-Hasselwander et al., 1979, Annals of the New York Academy of Sciences, Proc. of the 
Ninth Texas Symp., 226, 211 

Morselli A. et al., 2013, NuPhS, 239, 193 

Perotti E et al., 2006, NIM A, 556, 228 

Prest M., Barbiellini G., Bordignon G. et al., 2003, NIM A, 501, 280 
Swanenburg, B.N. et al., 1981, ApJ, 242, E69 
Tavani M., Barbiellini G., Argan A. et al., 2009, A&A, 502, 995 
Thompson D.J. et al., 1993, ApJS, 86, 629. 

Thompson D.J., et al., 1998, in Proc. 4th CGRO Symp., AlP Conf. Ser. n. 410, p. 39. 

Yang, H.J., Roe, B.P & Zhu, J., 2005, NIM A, 555, 370 


This preprint was prepared with the AAS macros v5.2. 


8. Appendix: The AGILE GRID vs. Fermi-LAT 


AGILE is an ASI Small Scientific Mission (Tavani et al. 2009) of total weight of 320 kg. It 
carries a scientific instrument dedicated to high-energy astrophysics whose heart is the Gamma- 
Ray Imaging Detector (GRID). Complementary items are the imaging Super-AGILE detector sen¬ 
sitive in the range 20-60 keV (Eeroci et al., 2007), the Mini-Calorimeter (Labanti et al. 2006), 
and the Anticoincidence system (Perotti et al. 2006). AGILE was launched in April 2007 and is 
operational in an equatorial orbit of average height of 530 km. 

Fermi is a NASA mission of a large class with a broad international collaboration. Its imaging 
gamma-ray instrument is based on a 16-unit Tracker (EAT) (Atwood et al., 2009), which is com¬ 
plemented by a massive Calorimeter and an Anticoindence system. Eermi was launched in June 
2008 in an orbit with inclination of 25° and average height of 550-600 km. 

Table Al summarizes the main parameters of the two instrument configurations, which are 
relevant for the angular resolution determination. 

Table Al: A comparison between the AGILE-GRID and Fermi-LAT 


Parameter 

AGILE-GRID 

Eermi-LAT 

Number of towers 

1 

16 

Total number of Tracker planes 

12 

18 

Vertical spacing (5) between adjacent planes 

1.8 cm 

3.2 cm 

Silicon tile size 

9.5 X 9.5 cm^ 

8.95 X 8.95 cm^ 

Silicon detector array for each plane 

4x4 

4x4 

Silicon-strip pitch (dp) 

121 yum 

228 yum 

Readout pitch 

242 yum 

228 yum 

Signal readout 

analog 

digital 

Ratio 6pl s 

0.007 

0.007 

Tungsten converter thickness per plane 

0.07 

0.03 Xo (front) 
0.18 (back) 

Number of planes with W converter 

10 

12 (front) 

4 (back) 

On-axis total radiation length 

0.9 

0.5 (front) 

0.8 (back) 

Total n. of readout channels 

36,864 

884,736 

Power consumption/channels 

400 yUW 

180yUW 





